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Method of Identifying Nucleic Acids 

Related Applications 

This application claims priority to USSN 60/1 15,109, filed January 8, 1999, which is 
incorporated herein in its entirety. 

Field of the Invention 

The present invention relates to nucleic acids and more particularly to methods of 
equalizing the representation of nucleic acids in a population of nucleic acid molecules. 

Background of the Invention 

Approximately 10,000-20,000 genes are thought to be expressed within living cells, 
depending upon the specific cell type. RNAs corresponding to different genes can be present in 
different levels in cells. For example, transcripts from as few as 10-15 genes may represent 10- 
1 5% of cellular mRNA by mass. In addition to these highly abundant transcripts, another 1 000- 
2000 genes encode moderately abundant transcripts, which can account for up to 50% of cellular 
mRNA mass. Transcripts from the remaining genes fall into the low abundance class. 

Because many genes are identified by isolating complementary DNA (cDNA) 
corresponding to an RNA sequence, a significant problem can arise because of differences in the 
levels at which specific RNAs are present in cell types. The most abundant sequences can be 
repeatedly sampled, while the lowest abundance class may be rarely, if ever, sampled. 

Several normalization and subtractive hybridization protocols have been developed to 
help overcome this problem. These techniques can be technically difficult to perform, and they 
can fail to detect cDNAs corresponding to rare transcripts. 

Summary of the Invention 

The invention is based in part on the discovery of novel procedures for equalizing, or 
normalizing, the representation of nucleic acids in a sample of nucleic acids in which different 
nucleic acids are initially present in the sample in unequal amounts. 
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Accordingly, in one aspect the invention provides a method of screening a population of 
nucleic acid sequences. The method includes providing a population of nucleic acid sequences, 
partitioning the population into one or more subpopulations of nucleic acids, and identifying a 
first nucleic acid sequence having an increased level in the subpopulation relative to its level in 
the starting population of nucleic acids. The first nucleic acid is then compared to a reference 
nucleic acid sequence or sequences. The absence of the first nucleic acid sequence in the 
reference nucleic acid or nucleic acid sequences indicates the first nucleic acid is a novel nucleic 
acid sequence. 

The RNA can be derived from a plant, a single-celled animal, a multi-cellular animal, a 
bacterium, a virus, a fungus, or a yeast. If desired, the RNA can also be partitioned prior to 
synthesizing cDNA. 

Among the advantages of the methods are that they eliminate, or minimize, redundant 
identification and characterization of identical nucleic acid sequences in a population of nucleic 
acids.. 

In some embodiments, the cDNA is synthesized to selectively generate cDNA species 
that are enriched for those sequences oriented towards the 5-terminus of the cDNA. In other 
embodiments, the cDNA is synthesized to enrich for those sequences oriented towards the 
3 ! -terminus of the cDNA. 

In some embodiments, the population is normalized by digesting the cDNAs with one or 
more restriction endonucleases, in different reaction vessels, so as to generate segregated 
multiple partitions. Preferably, each specific digested cDNA-fragment will occur in only one 
partition. 

In some embodiments, the cDNAs are partitioned by physical methods, which may 
optionally follow the restriction endonuclease digestion. The physical methods separate the 
cDNAs a function of their terminal nucleotide sequences, overall length and migratory pattern on 
a sizing matrix that possesses the ability to separate molecules as a function of their physical 
and/or biochemical properties. 
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In other embodiments, the cDNAs are partitioned during subsequent PCR-based 
amplification of adapter-ligated cDNA fragments that have been digested with one or more 
restriction endonucleases. 

In other embodiments, the cDNAs are partitioned by screening the original mixture of 
5 cDNAs so as to remove those sequences that have already been characterized. Screening occurs 
using partitioned subtraction, whereby the original cDNAs are brought into contact with a 
prepared, subtraction library of known sequence in such a way that any sequence contained 
within the original library that is complimentary to any element of the subtraction library is 
removed or suppressed. 

10 cDNA sequences may also be partitioned by determining the size of each cDNA fragment 

prior to sequencing; biasing for formation of larger fragment PCR products by lariat formation. 
In this method, a bias for the larger fragment within the PCR reaction is introduced to allow 
efficient preferential amplification of longer fragments. Alternatively, partitioning may occur by 
preferentially amplifying 5' terminal or 3' terminal sequences of mRNA molecules. 

15 If desired, the amplified cDNAs may fractioned by separating the amplified cDNAs on a 

sizing matrix that separates molecules as a function of their physical and/or biochemical 
properties and excising individual cDNA fragments from said sizing matrix. The excised cDNA 
fragments are then inserted into a recombinant vector, or further amplified. 

In some embodiments, the restriction endonuclease is a restriction endonuclease that 
20 possesses a recognition sequence 4 to 8 basepairs in length and produces either a 5 f - or 3-terminal 
overhang 0 to 6 basepairs in length. 

In some embodiments, the identified sequence is subjected to computational analysis. 
The computational analysis can include querying, or searching, a nucleotide sequence database to 
identify sequences that match, or the absence of any sequences that match. The database 
25 includes a plurality of known nucleotide sequences of nucleic acids that may be present in the 
sample. 
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Preferably, the nucleic acid database comprises substantially all the known, expressed 
nucleic acid sequences derived from a group comprising a plant, a single-celled animal, a multi- 
cellular animal, a bacterium, a virus, a fungus, or a yeast. 

In some embodiments, sizing includes diluting and re-amplification of the cDNAs, 
5 fractionating the re-amplified cDNAs by use of one or more sizing matrixes that separate the 
molecules as a function of their physical and/or biochemical characteristics, physically dividing 
or cutting the sizing matrixes into a plurality of sections, wherein each section is comprised of 
one or more cDNAs of similar molecular weight or size. The cDNAs are eluted from each of the 
sizing matrix section, ligated into a cloning vector and transformed into a host, e.g., a bacterial 
10 host. A plurality of the transformed host colonies are selected so as to ensure a statistically- 
accurate representation of the cDNAs originally contained within the sizing matrix sections. The 
□ inserts from this plurality of colonies are recovered and their molecular weight or size of are 
1 determined. A plurality of insert DNAs, wherein each successive insert has a molecular weight 
!^ or size that is within a 0.2 basepair window; and wherein only those DNA species that fall within 
W 15 the 0.2 basepair window is subsequently subjected to nucleotide sequencing. 

W As utilized herein, the term "normalized" is defined as a mixture of mRNAs (or cDNAs 

!=* thereof) in which the copy number of highly abundant mRNA species is reduced relative to its 
u copy number in a starting population of nucleic acids, and the copy number of a less abundant 
^ mRNA species has been enriched relative to the copy number of the latter mRNA in the starting 
'O 20 population. 

Among the advantages provided by the present invention are that it multiple partitioning 
strategies function in a synergistic manner so as to ameliorate unnecessary, redundant sequencing 
of the same sequence(s), while concomitantly enhancing the sequencing of rarer sequences. 

The partition strategies disclosed herein also normalize cDNA abundance by separating 
25 the cDNA sequences into multiple partitions possessing minimal sequence overlap. In addition, 
the various partitioning strategies are performed so as to assure that substantially all cDNAs are 
sampled. An additional normalization effect may be obtained by separating the resulting DNA 
fragments based upon their overall size (i.e., size fractionation). Moreover, it is also possible to 



4 



Attorney Ref.: 15966-534 



normalize the abundance of the cDNAs to an even greater degree by the use of one of several 
disclosed pre-characterization methods. 

All technical and scientific terms used herein have the same meanings commonly 
understood by one of ordinary skill in the art to which this invention belongs. Although any 
methods and materials similar or equivalent to those described herein can be used in the practice 
of the present invention, the preferred methods and materials are now described. The citation or 
identification of any reference within this application shall not be construed as an admission that 
such reference is available as prior art to the present invention. All publications mentioned 
herein are incorporated herein in their entirety by reference. 

Brief Description of the Drawings 

FIG. 1 is a flow diagram illustrating a method for normalizing the abundance of nucleic 
acid molecules in a population of nucleic acid molecules. 

FIG. 2 is a flow diagram illustrating a method of 5'-enriched cDNA synthesis according 
to the invention. 

FIG. 3 A is a schematic diagram showing restriction enzyme digestion and adapter 
ligation for enrichment of 5' ends of mRNA molecules. 

FIG. 3B is a histogram showing the regions of genes covered by clones constructed using 
5' end enrichment. 

FIG. 3C is a schematic diagram showing restriction enzyme digestion and adapter 
ligation for enrichment of mRNA molecules containing internal restriction fragments. 

FIG. 3D is a histogram showing the regions of genes covered by clones constructed 
using enrichment for internal restriction fragments. 

FIGS. 4 A and 4B are schematic illustrations showing the effects of partitioning on the 
types of nucleic acids recovered in relation to the abundance of the mRNA molecules. 
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Detailed Description of the Invention 

The present invention provides methods for identifying nucleic acids in a population of 
nucleic acid samples. It is based in part on normalizing the representation of sequences that may 
be initially present in different levels in the population of nucleic acid sequences. The 
5 normalization takes place by one or more methods of partitioning the nucleic acid population. 

A schematized overview of the invention is shown in FIG. 1 . At the input step 1 00 a 
starting population of RNA is chosen for analysis. Unless indicated otherwise, reference to a 
given RNA or population of RNAs is understood to also encompass reference to the 
corresponding cDNA or cDNAs. 
10 Any population of RNA molecules can be used as long as the population contains, or is 

suspected of containing, two or more distinct RNA molecules. The population can be isolated 

O i • i_ 

# from a starting sample using standard methods for isolating RNA. The RNA population can be 

u isolated from, e.g. , an entire organism or multiple organisms, or from a tissue or cell of an 

^ organisms. The RNA can also be isolated from, e.g. , cultured cells, such as eukaryotic or 

ffl 15 prokaryotic cells grown in vitro. If desired, the RNA can be mRNA, (e.g., polyA+ RNA), or 

J 8 stable RNAs (e.g. , ribosomal RNA, transfer RNA, or small nuclear RNA). The input RNA or 

2 cDNA can be a subpopulation containing the 5' end of RNA molecules (110), a subpopulation 

!^ having an internal regions of starting RNA molecules (1 12), or subpopulations containing the V 

yQ end ofthecDNA molecules (114). 

20 The selected population or subpopulation is next subjected to a normalization analysis 

(200). The normalization analysis includes one or more partitioning steps that decrease the 
relative amount of sequences that are abundant in the starting population of nucleic acids and 
increase the relative representation of sequences that are rare in the starting population of nucleic 
acids. A partitioning step can take place before or after mRNA is converted to cDNA. A 

25 partitioning step can also take place following amplification of a cDNA. Unless stated 

otherwise, any partitioning method described herein can be used in conjunction with one or more 
additional partitioning methods. Examples of suitable partitioning steps are provided below. 
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In some embodiments, cDNA molecules are subjected to digestion with restriction 
enzymes, after which adapter oligonucleotides are ligated to the digestion products, and the 
resulting products amplified. FIG. 1 indicates two types of digestions and adapter ligations 
which can be performed. The first, designated short chemistry (216) because it tends to result in 
5 shorter amplification products, uses two restriction enzymes, followed by ligation of adapter 
oligonucleotides having termini complementary to the termini of the internal digestion 
fragments. The second, designated long chemistry (218), similarly uses restriction digestion and 
adapter ligation but uses longer adapters, which generally result in longer amplification products. 

FIG. 1 also illustrates that the modified cDNAs can be subjected to size fractionation 
10 (220), which is an example of a partitioning method, and that information from the size fraction 
analysis can be used in a precharacterization analysis (222). A precharacterization can include, 
e.g., comparing the size of the insert to sequence databases of fragments sizes produced by the 
restriction enzyme. Amplification of short and long chemistry fragments can also be performed 
in association with partitioning steps, which are explained in detail below. 

15 The amplified products are next sequenced (300). Sequencing can be performed by any 

method known in the art. The compiled sequence data are then assembled (400), and the 
sequence generated is compared to known sequences, e.g., sequences in publicly available 
databases. 

The methods herein described are therefore useful for identifying genes, e.g., expressed 
20 genes in an organism of interest, e.g. , a human. The sequence information obtained is 

particularly useful for identifying genes transcribed at low levels, or generating low levels of 
steady state transcripts. The methods can also be used, e.g., to identify secreted proteins for 
potential therapeutic use and/or for drug targets; identify variations within the human genome, 
such as single nucleotide polymorphisms (SNPs); identify differences between normal and 
25 diseased tissue; and analyze differential gene expression in different tissues and/or species. 

Partitioning prior to cDNA synthesis 

One approach to normalize levels of mRNA from a given sample, e.g. a given cell or 
tissue type, is to arbitrarily separate a starting population of RNA molecules into many smaller 
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subpopulations, or collections. In general, a greater number of partitions increases the likelihood 
that a given partitions will lack a sequence or sequences that is abundant in the starting 
population of nucleic acid sequences. This method therefore allows for access to sequences that 
are expressed in very low copy number. 

Alternatively, RNA populations can be isolated from different cell types. This 
partitioning strategy is based on the premise that different tissues tend to express different 
subsets of genes. Thus, RNA sequences can be partitioned by sequencing multiple different 
cDNA libraries extracted from one or more tissues within the body. However, the partitioning 
will not typically be complete, because many genes are expressed in more than one tissue type. 

Synthesis and Amplification of cDNA molecules 

Typically, partitioning is performed on cDNA populations that have been modified for 
subsequent analysis. The modifications may include: (i) digesting the cDNA with at least one 
restriction endonuclease; (if) ligating an adapter oligonucleotide to one or more ends of the 
termini of the digestion products; and (Hi) amplifying the ligated products, e.g., in PCR-mediated 
amplification. These methods are particularly suited to cDNA molecule that have been 
constructed from the 5', internal, and 3' subpopulation of RNA molecules as described above. 
These manipulations are collectively known as SeqCalling™ chemistry. In preferred 
embodiments, cDNA is generated from populations of RNA molecules that have been divided 
into subpopulations containing 5' ends of transcripts, populations of molecules containing 
internal regions of RNA molecules, or subpopulations containing 3' ends of RNA molecules. 

A. Construction and amplification of cDNA subpopulation enriched for the 
5 ' ends of mRNA molecules 

5'-enriched cDNA synthesis generates cDNA species that are enriched for those 
sequences oriented towards the S'-terminus of the cDNA, and in which a specific oligonucleotide 
sequence is ligated to the 5 r -terminus. Approaches for generating cDNAs specifically enriched in 
transcript 5' ends are often based on the synthesis of a homopolymeric (e.g., dG or dA) tail by the 
enzyme terminal deoxynucleotidyl transferase (TdT) subsequent to the synthesis of the first 
cDNA strand. Second strand synthesis is then primed by the use of a complementary homo- 
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oligonucleotide primer sequence. See e.g., Frohman, et al, 1988. Proc. Natl. Acad. Sci. USA 85: 
8998-9002; Delort, et al, 1989. Nucl. Acids Res. 17: 6439-6448; Loh, et al, 1989. Science 243: 
217-220; Belyavsky, et al, 1989. Nucl. Acids Res. 17: 2919-2932; Ohara, et al, 1989. Proc. 
Natl. Acad. Sci. USA 86: 5673-5677. 

Alternatively, amplification can exploit the 5'-terminal cap structure present in eukaryotic 
mRNAs (see e.g., Furuichi & Miura, 1975. Nature 253: 374-375; Banerjee, 1980. Microbiol. 
Rev. 44: 175-205; Shatkin, 1985. Ce//40: 223-224). However, mRNA preparations generally 
include a mixture of both capped and non-capped mRNA species. The non-capped mRNAs are 
thought to be primarily the result of degradation within the cell or during the isolation procedure. 
An alternative approach to enrich for full-length mRNAs is to purify capped mRNA using 
affinity reagents. These reagents include naturally occurring proteins that bind the cap structure 
(see e.g., Edery, et al, 1995. Mol Cell. Biol. 15: 3363-3371); anti-cap antibodies (see e.g., 
Bochnig, et al, 1987. Eur J Biochem. 68: 460-467); and chemical modification of the cap, 
followed by selection for the modified cap structure (see e.g., Carninci, et al, 1996. Genomics 
37: 327-336). In addition, 5-oligo capping can also be used, in which specific oligonucleotide 
sequences are selectively added to 5'-capped mRNAs prior to first strand cDNA synthesis. 
Subsequent synthesis of the second strand, is primed by an oligonucleotide that is 
complementary to the modified cap sequence. See e.g., Maruyama & Sugano, 1994. Gene 138: 
171-174; Suzyki, et al, 1997. Gene 200: 149-156; Fromont-Racine, et al, 1993. Nucl. Acids Res. 
21: 1683-1684; U.S. Patent No. 5,597,713). 

An alternative method for isolating RNA molecules containing a capped 5' end is shown 
in FIG. 2. FIG. 2 depicts a flow diagram for 5'-enriched cDNA synthesis using a full-length 
mRNA having a 5 '-terminal cap sequence (Gppp) and a poly A+ tail. Also shown in FIG. 2 is 
truncated mRNA having a 5' terminal phosphate group. Typically, RNA preparations contain a 
mixture of full-length capped RNAs and truncated mRNAs. The truncated RNAs can arise, e.g. , 
by intracellular degradation of the RNA or by degradation of the RNA during its isolation. 

In the first step in FIG. 2, the free 5'-terminal phosphate groups of the truncated or 
degraded mRNAs are removed by the action of a phosphatase, e.g., the bacterial alkaline 
phosphatase shown, or calf intestinal alkaline phosphatase. The phosphatase is then inactivated. 
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In the second step, the 5' cap is removed from the full-length mRNA using a pyrophosphatase, 
e.g., the tobacco acid pyrophosphatase shown in FIG. 2. The resulting product is the decapped 
full-length RNA with a free 5'-terminal phosphate group. 

In the third step in FIG. 2 , the phosphate group serves as a substrate for an RNA ligase- 
5 mediated reaction that attaches a specific DNA/RNA hybrid to the 5'-terminus of the full-length 
mRNAs. An RNA containing the ligated hybrid is used as a substrate for first and second strand 
cDNA synthesis. Preferably, a combination of oligo(dT)- and random hexamer-mediated first 
strand priming is performed in the presence of E. coli ligase to enhance overall cDNA length. 
Preferably, an RNase and thermal cycling are used to remove the RNA strand after first strand 
10 synthesis. The resulting single strand DNA (ssDNA) functions as a more effective reagent for 
the priming of second strand synthesis. 

Although first strand synthesis occurs for both types of mRNA species {i.e., full-length 
and truncated/degraded), only those mRNAs with the appropriate sequence ligated to the 5'- 
terminus (i.e., full-length mRNAs) contain a priming site for subsequent second strand synthesis. 
1 5 Thus, RNAs derived from the full-length mRNAs are selectively amplified. 

Preferably, a thermostable enzyme for second strand synthesis in a non-thermal cycled 
temperature profile is used to ensure more stringent priming of the second strand reaction 
compared to a non-thermostable enzyme. 

A double-stranded cDNA prepared with an adapter containing an oligonucleotide 
20 sequence (nR plus "signature sequence") ligated to the 5'-terminus is digested with a restriction 
endonuclease as shown in FIG. 3 A. The oligonucleotide RS [SEQ ID NO:l] (or nR) is used to 
prime the PCR amplification step subsequent to the ligation of the restriction digestion products. 
The nJ/nJ PCR product is shown as lined-through to denote that it does not clone efficiently in E. 
coli. 

25 A representation of the distribution of clones derived using 5' enriched synthesis with 

respect to the region of the gene they include is shown in FIG. 3B. A reference mRNA 
containing a 5' terminus, an ATG initiation codon, a Stop codon, and a 3' terminus is shown 
along the X-axis. Also shown is a histogram showing the number of clones (Y-axis) containing 
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sequences derived from the indicated regions of the reference mRNA. The histogram reveals that 
the 5' enrichment method method generates distributions enriched in 5' end fragments, and has 
increased proportions of fragments containing the start codon and the adjacent 90 bp of coding 
sequences. 

B. Construction and amplification of cDNA subpopulations enriched for the 
interior regions ends of RNA molecules 

To generate relatively short cDNA fragments generated from the interior regions of a 
RNA molecule, i.e., from a region not containing the 5' or 3 ' terminus, the following procedure 
is used. 

RNA is purified using any standard procedure (see e.g., Berger, 1987. Methods Enzymol 
152 : 215-219) and cDNA is synthesized according to standard protocols, such as random 
oligomer or oligo-dT primed synthesis (see, e.g., Gubler & Hoffman, 1983, Gene 25: 263-269, 
Okayama& Berg, 1982, Mol. Cell Biol. 2: 161-170). 

The cDNA is initially digested with a pair of restriction endonucleases. Although any 
enzyme pair that generates distinct 5 T -terminus overhangs is acceptable, a preferred embodiment 
utilizes enzymes that possess a 4-8 basepair (bp) recognition site yielding a 0-6 bp S'-terminal 
overhang, and a more preferred embodiment utilizes enzymes that possess a 6 bp recognition 
sequence and generates a 4 bp 5 f -terminus overhang. One form of manipulation for generating 
internal fragments is shown in FIG. 3C. The cDNAs are digested with two restriction 
endonucleases, yielding three types of fragments (two "homo", one "hetero" termini). Following 
digestion, specific adapters are ligated and the fragments are PCR amplified based upon the 
specific adapter sequence utilized. As indicated by the crossed lines, the nR-nR and nJ-nJ 
fragments are unstable in E. coli, and are rarely observed following cloning. 

Two suitable 24 nucleotide adapter molecules can be generated from RA24 [SEQ ID 
NO:9]; RC24 [SEQ ID NO:10]; JA24 [SEQ ID NO:l 1]; or JC24 [SEQ ID NO: 12]. The adapters 
are generated by annealing the RA24, RC24, JA24 or JC24 24-mer oligonucleotides [SEQ ID 
NOs:9-12, respectively] with 12-mer oligonucleotides possessing sequences that are 
complementary to the last 8 nt of the 3'-terminus of the 24-mer and the 4 bp overhang. The 
sequences of these primers and other primers described herein are provided in Table 1. 
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These 4 bp overhang sequences are chosen so as to be complementary to the overhangs 
that are generated by the restriction endonuclease digestions. In addition, the last 3'-terminal 
nucleotide of the 24-mer adapter (i.e., A or C) is selected such that a functional restriction 
endonuclease recognition site is not re-generated when the adapter anneals to the digested cDNA. 

5 Following ligation of the adapters, the restriction endonucleases are heat-inactivated, and 

the reaction mixture is PCR amplified. 

Internal fragments may alternatively be generated using a second type of adapters, which 
results in longer amplified fragments (also referred to as "Long Internal Chemistry" or "Long 
Chemistry"). This method is similar to short chemistry, except all adapters possess an additional 
10 common sequence on their 5'-termini. This technique suppresses the amplification of small 

fragments while concomitantly increasing the amplification of longer fragments. The subsequent 
5 PCR amplification with the "X" and "J" primers results in production of both a hetero (/. e. , "RX- 
;* -JR") adapter fragment and "homo" adapter fragments (/. e., "RX-XR" and "RJ-JR"), which are 
y unstable in a host and are rarely observed following the cloning process. 

m 15 The effectiveness of enriching for internal fragments is shown in FIG. 3D. Several 

I thousand sequences generated from internal cDNA framents and compared against a database of 

:~ approximately 5000 known genes with annotated start and stop sites. Each sequence matching 
^ the database was assigned a location on the gene relative to the start (0.0) and stop (1 .0) locations 
fi relative to the location of the 5'-most matching nucleotide (of the gene). The distribution from a 
^ 20 standard run shows that most fragments are located "internally" (i.e., within the coding region). 
Fragments covering the start codon plus an additional 90 bp (located immediately 3' of the start 
codon) are significant, because they have a high probability of containing enough sequence to 
identify secreted proteins. A small but significant fraction of the fragments covers the start 
codon and the additional 90 bp. 

25 Following digestion, adapters are ligated to these 5'-terminal overhangs. The primers are 

longer relative to primers used to generate short fragments. Two specific pairs of adapter 
molecules that can be used in long chemistry synthesis include RXC [SEQ ID NO:2]; RXA 
[SEQ ID NO:3]; RJC [SEQ ID NO:4]; or RJA [SEQ ID NO:5]. The adapters are generated by 
annealing RXC, RXA, RJC or RJA oligonucleotides [SEQ ID NOs:2-5, respectively] with 12- 
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mer oligonucleotides possessing sequences that are complementary to the last 8 nt of the 
3*-terminus of the 24-mer and the 4 bp overhang. These 4 bp overhang sequences are chosen so 
as to be complementary to the overhangs that are generated by the restriction endonuclease 
digestions. In addition, the last 3-terminal nucleotide of the 24-mer adapter (i.e., A or C) is 
selected such that a functional restriction endonuclease recognition site is not re-generated when 
the adapter anneals to the digested cDNA. 

Following the ligation of the adapters, the restriction endonucleases are heat inactivated 
and the reaction mixture is PCR amplified. While the sequences of the two adapters are distinct, 
they nevertheless possess common 5' sequences that allow the formation of lariat or pan-handle 
structures that function to suppress PCR-mediated amplification of the shorter fragments. 

C. cDNA Synthesis of molecules enriched for 3 ' ends 

3'-enriched cDNA synthesis generates cDNAs that are enriched for the sequences 
oriented towards the 3'-terminus of the cDNA. This is accomplished by synthesis of the first- 
strand using a specific oligonucleotide sequence that has been modified to contain an adapter 
sequence at its 5-terminus [SEQ ID NO: 14]. Following first-stand cDNA synthesis with the 
primer, standard cDNA synthesis protocols are utilized as illustrated in FIG. 2. 

The 3-enriched cDNA is digested with one restriction endonuclease. Although any 
enzyme that generates a distinct 5'-terminus overhang is acceptable, it is generally most preferred 
to utilize an enzyme that possesses a 6 bp recognition site yielding a 4 bp 5'-terminal overhang. 
Following digestion, an adapter is then ligated to these 5'-terminal overhangs. These adapters are 
generated from the JA24 [SEQ ID NO:l 1] or JC24 [SEQ ID NO: 12] 24-mer annealed with 
12-mer oligonucleotides possessing sequences that are complementary to the last 8 nt of the 
3-terminus of the 24-mer and the 4 bp overhang. These 4 bp overhang sequences are chosen so 
as to be complementary to the overhangs that are generated by the restriction endonuclease 
digestions. In addition, the last 3'-terminal nucleotide of the 24-mer adapter (i.e., A or C) is 
selected such that a functional restriction endonuclease recognition site is not re-generated when 
the adapter anneals to the digested cDNA. 
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Following the ligation of the adapters, the restriction endonucleases are heat inactivated 
and the reaction mixture is PCR amplified 

Longer fragments enriched for the 3'-ends can be obtained by ligating a longer primer to 
cDNA molecules that have been digested with a restriction enzyme. Any enzyme that generates 

5 a distinct 5 r -terminus overhang can be used. It is generally preferred to utilize an enzyme that 
possesses a 6 bp recognition site yielding a 4 bp S'-terminal overhang. Following digestion, an 
adapter is then ligated to the 5 -terminal overhangs. Acceptable adapters are generated from the 
JA24 [SEQ ID NO:l 1] or JC24 [SEQ ID NO: 12] 24-mer annealed with 12-mer oligonucleotides 
possessing sequences that are complementary to the last 8 nt of the 3'-terminus of the 24-mer and 

10 the 4 bp overhang. These 4 bp overhang sequences are chosen so as to be complementary to the 
overhangs that are generated by the restriction endonuclease digestion. In addition, the last 3- 
terminal nucleotide of the 24-mer adapter (i.e., A or C) is selected such that a functional 
restriction endonuclease recognition site is not regenerated when the adapter anneals to the 
digested cDNA. 

15 While the sequences of the two adapters are distinct, they possess common 5 f sequences 

that allow the formation of structures that suppress PCR-mediated amplification of the shorter 
fragments. 

Following the ligation of the adapters, the restriction endonucleases are heat inactivated 
and the reaction mixture is PCR amplified. 

20 The cDNA fragments prepared as above can be size-fractionated, e.g. , electrophoretic 

fractionation on agarose or polyacrylamide gels, or other types of gels comprised of a similar 
material. The cDNA fragments may then be physically excised in defined size ranges (i.e., as 
identified by size makers) and recovered from the excised gel fragments. Additionally, if the 
quantities of isolated cDNA fragments are low, they can be amplified, e.g., by PCR amplification 

25 For example, if the cDNA fragments are generated by Long Internal SeqCalling™ Chemistry 
protocol, they are amplified with J23 [SEQ ID NO:6] and X22 [SEQ ID NO: 15] primers (either 
before or after fractionation) prior to cloning, as these cDNAs cannot be efficiently cloned into E. 
coll Similarly, if the cDNA fragments are generated by Long 5 ! SeqCalling™ Chemistry 
protocol, they can be amplified by J23 [SEQ ID NO:6] and RS [SEQ ID NO: 1] oligonucleotides 
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(either before or after fractionation) prior to cloning, as these products cannot be efficiently 
cloned into E. coli. 

When PCR amplification is used to amplify fragments, conditions are preferentially 
chosen to minimize non-productive hybridization events. It has been observed that DNA re- 
5 hybridization during the PCR amplification process (designated the "Cot effect"; see e.g., 

Mathieu-Daude, et al, 1996. Nucl Acids Res. 24: 2080-2084) can inhibit amplification. This 
effect is particularly evident during later PCR amplification cycles, when a substantial 
concentration of the amplified product has accumulated and the primer concentration has been 
depleted. As a result, amplification in the later PCR cycles typically follow non-linear dynamics. 

10 By manipulating PCR amplification reaction conditions, it is possible to markedly 

enhance the "Cot effect 1 ', by the insertion of a slow-annealing step in between the denaturation 
and re-naturation steps in each PCR amplification cycle. The slow-annealing temperature is 
chosen so as to be above that of the primer-template melting temperature (T m ), but at or above 
that of the template-template T m , thus favoring template-template annealing over template-primer 

15 annealing. For example, a 85-75°C decrease in temperature at a 10°C/minute gradient can be 
utilized 



Partitioning methods 

One or more of the following techniques, or combinations these techniques, can be used 
20 to normalize the abundance of RNA (or their cDNA counterpart) species within a given cell or 
tissue sample. 

(i) Partitioning by restriction endonuclease digestion 

A cDNA library can be partitioned into many different sets of fragments by digestion 
with different restriction enzyme pairs. Fragmentation of the same cDNA library with different 
25 sets of restriction enzymes, in different reaction vessels, results in segregated multiple partitions, 
i.e., each specific fragment will occur in only one partition. The digested fragments can be 
analyzed further, e.g., by direct sequencing, cloning of the digested fragments or sequencing, or 
one or more of these techniques. 
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If desired, the cDNA is digested into fragments of a length that is convenient for 
sequencing. Preferably, multiple different partitions, e.g., 10-100, 20-750, or 50-250 partitions 
are obtained. 

5 (H) Partitioning by fragment size or other physical property 

Partitioning can also be performed using other separation methods that separate DNA 
molecules according to their physical characteristics. The methods can include, e.g., separation 
based on physical and/or biochemical properties {i.e., molecular weight/size, terminal nucleotide 
sequences, exact migratory pattern, and the like). Separation methods can include, e.g., gel 
10 electrophoresis, including agarose or polyacrylamide gel electrophoresis, high pressure liquid 

chromatography (HPLC), preparative-scale capillary electrophoresis, and similar methodologies. 

K U In one embodiment, unique cDNAs that represent unique (i.e., not previously sequenced) 

U fragments are selected based on their presence in a characteristic restriction enzyme fragment. In 

; ~1 this process, a cDNA population is digested with restriction endonucleases, fractionated, and 

! J5 15 fragments in a desired size range are recovered. The recovered fragments are then ligated to a 

* vector and transformed into an appropriate host, e.g., E. colu Rather that being directly 

S sequenced following the selection process, the DNA fragments are isolated and separated, e.g., 

t*\ sized using one or more sizing matrixes that separate the molecules as a function of their physical 

0 or biochemical properties. The embodiment is thus referrred to as "clone sizing". Those 

^ 20 recombinant clones that have an insert with characteristics not present in a reference database are 

determined to contain a unique DNA fragment. Preferably, only unique fragments are 

subsequently sequenced. 

For example, a DNA fragment that is sized in this way possesses two pieces of 
information that serve as a unique identifier: (/) the identity of the restriction endonuclease used 
25 to generate the fragment, and (//) the size of the fragment. With these two pieces of information, 
fragments are picked for subsequent nucleotide sequencing by searching for a specific fragment 
within a 0.2 basepair window. If a fragment is present in the window, the E. coli clone 
containing the fragment is re-arrayed on a liquid handling robot such as a Tecan Genesis or 
Packard Multiprobe device, and sequenced. When multiple fragments are present within the 0.2 
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bp window, only one is selected to be sequenced. Thus, by use of this sizing filter, sequencing of 
identical fragments is significantly lowered. 

By sizing individual fragments and comparing the observed size to previously determined 
sequences, i.e., using a "sizing filter", only fragments of unique lengths need to be sequenced. 

5 To pre-size large numbers of fragments, the fragments can be initially pooled as a 

function of their expected size, so as to ensure the any fragment occurs in a minimum of at least 
three individual pools. 

Size fractionation may be accomplished in a number of ways. One commonly utilized 
method is electrophoretic fractionation on agarose or poly aery lamide gels, or other types of gels 
10 comprised of a similar material The cDNA fragments may then be physically excised in defined 
size ranges (i.e., as identified by size makers) and recovered from the excised gel fragments. 
J Additionally, if the quantities of isolated cDNA fragments are low, they can be PCR amplified at 
* this stage. For example, if the cDNA fragments are generated by Long Internal SeqCalling™ 
H Chemistry protocol, described above, they must be amplified with J23 and X22 primers (either 
m'15 before or after fractionation) prior to cloning, as these cDNAs cannot be efficiently cloned into E. 
|? 5 coll Similarly, if the cDNA fragments are generated by Long 5' SeqCalling™ Chemistry 
f protocol, described above, they must be amplified by J23 and RS oligonucleotides (either before 
or after fractionation) prior to cloning, as these products cannot be efficiently cloned into E. coll. 
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(Hi) Partitioning based on hybridization 

Screening can be performed using a variety of methods that rely on hybridization 
between a probe sequence or sequences and a cDNA library. Members of the library containing 
a homologous sequence are then removed from the library. For example, a cDNA library can be 
brought into contact with a prepared library of known sequence in such a way that any sequence 
contained within the substrate library that is complimentary to any element of the subtraction 
library is removed or suppressed. This method obviates re-characterizing, e.g., re-sequencing, 
already characterized members of the cDNA population. 

(iv) Amplification-associated partitioning 

Partitioning can also be performed in association with amplification. In particular, 
partitioning can be carried out during PCR amplification of adapter-ligated cDNA fragments 
described above. During PCR-mediated amplification of mixtures of cDNA fragments, short 
fragments tend to be preferentially amplified relative to large fragments. PCR conditions can be 
adjusted to favor the formation of larger fragments within the PCR reaction to allow efficient 
preferential amplification of longer fragments. 

Normally, two different primers are used in PCR amplification to prime the enzymatic 
activity of the polymerase at each terminus of the target sequence. Conversely, if primers with 
identical 5' sequences are used, there is a tendency for the fragments to form lariat or pan-handle 
structures, due to intra-strand hybridization, which interferes with the amplification process. 
Because the probability of the two ends of a polymer (i.e., cDNA fragment) finding one another 
is inversely proportional to a fractional power of the polymer length, short fragments tend to 
form these lariat structures more readily than do longer ones. Accordingly, this effect is 
exploited in the amplification of long cDNA fragments. See U.S. Patent No. 5,565,340, whose 
disclosure is incorporated herein by reference, in its entirety. 

Long fragment amplification can be enhanced using DNA fragments to which have been 
ligated long adapter sequences as described above. Amplification is dependent upon a number of 
factors that can alter the ratio of a linear adapter structure, which is permissive for amplification, 
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and a lariat-loop structure, which suppresses amplifications. The equilibrium constant associated 
with the formation of the suppressive and the permissive structures, and, therefore, the efficiency 
of suppression of particular DNA fragments during PCR, is primarily a function of the following 
factors: (/) differences in melting temperature of suppressive and permissive structures; (zi) 
5 position of the primer sequence within the adapter; (Hi) the length of the target DNA fragments; 
(iv) PCR primer concentration; and (v) primary structure. 

Analysis of partitioned cDNA molecules 

Partitioned cDNA molecules are next analyzed by comparing the sequences to a reference 
nucleic acid or nucleic acids. To facilitate analysis of partitioned cDNA molecules, they can, if 
10 not subcloned previously, be ligated into an appropriate vector and transformed into cells by any 
applicable method. 

The reference nucleic acid or nucleic acids can be any fragment for which sufficient 
T information is available to unambiguously identify the partitioned cDNA molecule. The 
H reference nucleic acid or nucleic acids can therefore be part of, e.g. , sequence databases, or 
S3 15 - databases of other characteristics that unambiguously identify a nucleic acid. Examples of such 

characteristics include e.g. , a compilation of fragment sizes associated with specific restriction 
^ enzymes for a particular gene. In some embodiments, partitioned nucleic acids will be 

sequenced. The partitioned sequences can be sequenced by any method known to the art and the 
2 resulting sequence data is analyzed by computer-based systems. 

20 Suitable databases include publicly available databases that comprehensively record all 

observed DNA sequences. Such databases include, e.g., GenBank from the National Center for 
Biotechnology Information (Bethesda, Md.), the EMBL Data Library at the European 
Bioinformatics Institute (Hinxton Hall, UK) and databases from the National Center for Genome 
Research (Santa Fe, N.Mex.). However, any database containing entries for the sequences likely 

25 to be present in such a sample to be analyzed is usable in the further steps of the computer 
methods. Methods of searching databases are described in detail in e.g., U.S. Patent No. 
5,871,697, whose disclosure is incorporated herein by reference, in its entirety. 
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Table 1 below summarizes the various primers and adapters disclosed herein. 



Table 1 



SEQ 


Name 


Sequence (from 5' to 3') 


1 
1 


rvo 


CTCTCCGATG 


CAGGTGGC 








z 


IvAL 


AGCACACTCC 


AGCCTCTCTC 


CGAGCACATG 


CGACACTGAG 


TACT AC 

J_ -LA. \ * _L_ ±± \ - 






AGCACACTCC 


AGCCTCTCTC 


CGAGCACATG 


CGACACTGAG 


TACTAA 


A 




AGCACACTCC 


AGCCTCTCTC 


CGAACCGACG 


TCGAATATCC 


ATGCAGC 






AGCACACTCC 


AGCCTCTCTC 


CGAACCGACG 


TCGAATATCC 


ATGCAGA 


6 


J23 


ACCGACGTCG 


AATATCCATG 


CAG 






7 


R23 


AGCACACTCC 


AGCCTCTCTC 


CGA 






8 


NR17 


AGCACACTCC 


AGCCTCT 








9 


RA24 


AGCACACTCC 


AGCCTCTCTC 


CGAA 






10 


RC24 


AGCACACTCC 


AGCCTCTCTC 


CGAC 






11 


JA24 


ACCGACGTCG 


AATATCCATG 


CAGA 






12 


JC24 


ACCGACGTCG 


AATATCCATG 


CAGC 






13 


dT-R 


AGCACACTCC 


AGCCTCTCTC 


CGA 






14 




AGCACACTCC 


AGCCTCTCTC 


CGATTTTTTT 


T T T T T T T T T T 


TTT 



EXAMPLES 

The invention will be further described in the following examples, which do not limit the 
scope of the invention described in the claims. Examples 1-6 collectively describe the synthesis 
and amplification of cDNA subtractions enriched for the 5 ? terminal sequences of mRNA 
molecules. Example 7 describes clone sizing. 

Example 1. 5 f cDNA Synthesis — phosphatase/pyrophosphate digestion 

For each reaction, 2.5 |ig mRNA (do not exceed 3 jag total) is added to H 2 0 so as to 
provide a total volume of 73.5 |iil. This mixture is then heated to 65°C for 10 minutes, and 
quick-cooled on ice. The CIAP Cocktail (see below) is made as follows: 

CIAP Cocktail : 

For each reaction: 1 0 jjlI 1 Ox CIAP buffer 1 1 0 (il 

2.5 |il RNasin (Promega) x 1 1 27.5 |il 

10^1 0.1 MDTT 110 *il 
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4^1 0.01 U/idCIAP* 35 pi 

1) 26.5 |J of the above enzyme mixture is added to each 3 \x\ mRNA to give 
a total volume of 30.5 (il. 73.5 jal of the RNA mix is then added to give a 



5 final volume of 1 00 jil. 





2) 


Incubate at 37°C for 40 minutes. 




3) 


Add 100 ul TE buffer (10 mM Tris pH 8.0; 0.1 mM EDTA). 




4) 


Add 200 ul Acid-Phenol. 




5) 


Mix vigorously. 


10 


6) 


Add 200 ul Chloroform-Isoamyl Alcohol (24:1 v/v). 




7) 


Mix vigorously. 




8) 


Centrifuge in a microfuge at maximum speed for 10 minutes. 




9) 


Remove supernatant and transfer to new tube. Discard bottom layer. 




10) 


Repeat steps 4-9 (only for CIAP treatment, not in later steps). 


a i5 


11) 


Add 2 ul ssDNA carrier and 20 ul 3 M Sodium Acetate to each tube. 




12) 


Vortex 10 seconds and add 440 ul of absolute ethanol. 




13) 


Vortex 10 seconds and incubate at least 30 minutes at -80°C. 




14) 


Centrifuge samples at 13,200 x g for 15 minutes. 




15) 


Wash nucleic acid pellets with 70% ethanol and air-dry pellet. 


20 


16) 


Dissolve nucleic acid pellet in 70 ul water and cool on ice. 




17) 


Centrifuge for 10-15 seconds at maximum speed. 




18) 


Transfer contents of tubes to 8-strip tubes. 




19) 


Add 30 ul TAP cocktail (see below). 
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TAP Cocktail : 

For each reaction: 1 0 pi 1 Ox TAP buffer 110 

2.5plRNasinxll 27.5 p.1 

15.5 pi H 2 0 170.5 pi 

2.0 jal 1 0 U/pl TAP (Epicenter) 22 pi 



20) Add 30 pi of above mixture to each 70p.l CIAP-treated sample for 
a total volume of 100 pi. 

2 1 ) Incubate at 3 7°C for 45 minutes. 

10 22) Repeat Phenol/Chloroform extraction and precipitation as above in 

steps 6-9 and then 1 1 - 1 5 (do not resuspend pellet). 
P Example 2. 5' cDNA Synthesis: DNA-RNA Hybrid Primer Ligation 

1) Transfer samples from Example 1 to 8-strip tubes. 
; ~1 2) Resuspend pellet in Ligation Cocktail (see below). 

3 15 

, 3 Ligation Cocktail : 

M For each reaction: 3 pi 10 mM ATP 33 pi 

Q lplRNasinxll 11 pi 

^ 4.5 pi FI 2 0 49.5 pi 

S 20 2 pi R-B AP-TAP DNA/RNA hybrid oligomer 22 pi 



3) Add 10.5 pi of above mixture to each pellet, dissolve pellet completely at 
room temperature by (preferably) tapping the tube or vortexing if needed. 

4) Make an enzyme mix as follows: 

25 Enzyme Mixture : 

For each reaction: 30 pi H 2 0 3 30 pi 

1 2 pi 5x DN A Ligase Buffer (Life Tech) x 1 1 1 32 pi 

1.5plRNasin 16.5 pi 

6 pi T 4 RNA Ligase (Life Tech.) 66 pi 

30 
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Total reaction volume 60 pi 

5) Incubate overnight at 20°C. 

6) Repeat Phenol/Chloroform and precipitation as above in CIP/TAP Cocktail 
protocol steps 6-9 and 11-15 (do not resuspend pellet). 

Example 3, 5 f cDNA Synthesis: cDNA First-Strand Synthesis 

1) Resuspend cDNA pellet in Random Hexamer Cocktail (see below). 

Random Hexamer Cocktail : 

For each reaction: 1 0 pi H 2 0 x 1 1 1 1 0 pi 

0.5 pi random hexamer (dN 6 -5'-Phosphate, 100 pM) 5.5 pi 
5 |il Oligo-(dT) (dT 30 VN-5Thosphate, 100 pM) 55 pi 

2) Add 15.5 pi of above mixture to each tube and resuspend pellet. 

3) Heat at 70°C for 10 minutes and quick-cool on ice. 

4) Make First-Strand Synthesis Cocktail as follows (see below). 

First-Strand Synthesis Cocktail : 



For each reaction: 6 pi 5x First-Strand Buffer 66 pi 

3pll0mMdNTPs 33 jal 

3pll00mMDTTxll 33 pi 

1 |nl RNase Inhibitor 1 1 pi 



5) Add 13 pi of the above mixture to each 15.5 pi sample to give a total volume 
of 28.5 pi. 

6) Incubate at 37°C for 2 minutes. 

7) Add 1 .5 pi Superscript II RT to each reaction for a total volume of 30 pi. 

8) Incubate at 37°C for 1 0 minutes. 

9) Incubate at 42°C for 1 hour. 
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10) Incubate at 16°C. 

1 1) Add 40 |iil of the following DNA Ligase Mixture (see below) to each reaction 
tube for a total volume of 70 jal. 

E. coli DNA Ligase Mixture : 

For each reaction: 4 (il 1 Ox E. coli Ligase Buffer x 1 1 44 \il 

33 jLtl H 2 0 330 (al 

3 p.1 E coli DNA Ligase (10 U/^l) 33 nl 

12) Continue incubation at 16°C for 2 hours. 

Example 4. 5 f cDNA Synthesis: removal of non-ligated Primers 

While the above 2 hour incubation described in Example 3 is progressing, prepare 
one Boehringer-Mannheim Quick-Spin G-50 columns per reaction as follows: 

1) Mix the resin bed well by inverting the columns repeatedly. 

2) Remove the top cap first, and then the bottom cap. This avoids bubble formation 
and resultant poor performance of the spin-column. 

3) Stand column vertically and allow to drain completely. 

4) Add 0.75 ml of 10 mM Tris (pH 7.5) to the top of the bed without disturbing. 
If the bed becomes disturbed, pipette the solution up and down slowly to mix 
the bed uniformly and allow the bed to re-settle so as to form a uniform surface. 

5) Stand column vertically and allow to drain completely. 

6) Place the columns into a 1 5 ml conical centrifuge tube with the vendor's 
associated collector tube beneath the spin-column to collect the sample. 

7) Centrifuge spin-column at 1000-1200 x g for 2 minutes. 

8) Remove spin-column with a forceps and remove the tube with flow through 
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and discard. 

9) Carefully load the sample to the top center of the spin-column. 

10) Wash the sample tube with 20 pi H 2 0 and load on the same column. 

1 1) Place a new collection tube beneath each spin-column and centrifuge at 
1000-1200 x g for 4 minutes. 

12) Remove spin-columns and collect the flow-through into new, labeled tubes. 

13) Total sample volume will be approximately 105 pi. 
Example 5. 5' cDNA Synthesis: RNase (H, A, and T t ) Treatment 

1) To each reaction described in Example 4 add Second-Strand Reaction Buffer (see 

below). 

Second-Strand Reaction Buffer : 

For each reaction: 3 pi 1 00 mM DTT 33 pi 

6 pi First-Strand Buffer 33 pi 

30 pi Second-Strand Buffer x 1 1 330 pi 

6 pi H 2 0 66 pi 



2) Add 45 pi of the above mixture to each 105 pi sample to give a total volume of 
150 pL 

3) Add 2 pi of RNase H to each sample. 

4) Incubate at 37°C for 30 minutes to nick the RNA in RNA/DNA hybrids. 

5) Make an RNase Mixture comprising: 22pl RNase H, 44 pi RNase Cocktail 
(Ambion; available as an RNase A and RNase T 1 mixture). 

6) Heat samples to 95°C for 2 minutes. 

7) Slow cool down to 37°C and continue incubation. 

8) Add 3 pi RNase Mixture to each of the cDNAs, mix by pipetting up and down. 
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9) Continue incubation at 37°C for an additional 10 minutes. 

10) Heat samples to 95 °C for 2 minutes. 

1 1) Slow cool down to 37°C and continue incubation. 

12) Add an additional 3 pi of RNase Mixture to each of the cDNAs, mix by 
5 pipetting up and down. 

13) Continue incubation at 37°C for an additional 15 minutes. 

14) Repeat Phenol/Chloroform extraction and precipitation as above in 
steps 6-9 and then 11-15. 

1 5) Dissolve pellet in 20 pi H 2 0. 

10 16) Remove a 5 pi aliquot for Second-Strand (see below) synthesis for producing 

5'-cDNA for SeqCalling™ Chemistry Protocol. 

Example 6. Second-Strand Synthesis for Producing S'-cDNA for SeqCalling™ Chemistry 

1) Generate PCR Mixture (see below) as follows: 

PCR Mixture: 

1 5 For each reaction: 5 pi 1 Ox PCR Buffer x 1 1 55 pi 

1 jliI 1 0 mM dNTPs 5.5 pi 

1 pi 10 jaM R17 Primer 5.5 jil 

37.5jalH 2 0 412.5 jal 

0.5 pi Advantage Polymerase 5.5 pi 

20 



2) Add 45 pi of the above mixture to each 5 pi sample, for a total volume 50 pi. 

3) Heat samples as per protocol below, making sure that the sample tubes are placed 
in the thermocycler only after it has reached >80°C. 

94°C for 2 minutes | 
25 55°C for 2 minutes | x 1 Cycle ONLY 

72°C for 60 minutes | (Cycle designated KM-AD-2N) 
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4°C for long-term storage 

4) Warm reaction tubes to 37°C. 

5) Make SAP Cocktail (see below) as follows 

SAP Cocktail : 

For each reaction: 1 2 jjl 1 Ox SAP Buffer x 1 1 1 32 (il 

5 nl H 2 0 55 jal 

3 |il Shrimp Alkaline Phosphatase (SAP; 1 U/jil) 33 pi 



6) Add 20 [il of SAP Cocktail to each reaction. 
10 7) Heat to 37°C for 30 minutes. 

8) Purify samples by Qiagen 96-well plate as manufacture's protocol. 
D 9) Elute cDNAs in 100 jlxI lOmM Tris-HCl buffer and proceed with fluorometry. 

:^ Example 7. Clone Sizing 

S SeqCalling™ Chemistry products generated in any of Examples 1-6 are diluted and re- 

y "45 amplified. Fractionation is then performed by electrophoresising the re-amplified sample on an 
M» agarose gel using MetaPhor agarose (FMC). After the electrophoresis, the gel is physically cut 
~T into a total of 48 fractions. 24 of the fractions are derived from a 4% MetaPhor gel, and 

correspond to the lower molecular weight fractions; whereas the other 24 fractions derived from 
J3 the 3% MetaPhor gel, correspond to the upper molecular weight fractions. 

20 Following the elution of the DNA from the gel fractions, the DNA fragments are ligated 

into a vector with the TOPO-TA cloning vector (Invitrogen). These plasmids are then 
transformed into E. coli. The transformed bacterial cells are plated onto petri dishes and grown 
to a size that allows automated colony picking. A suitable number of colonies/fraction are 
selected so as to ensure a statistically accurate representation of the DNA fragments contained 

25 within the fraction {i.e., suitable numbers of picked colonies/fraction are 48 or 96). Following 
the incubation of the selected clones, the fragment contained within each individual clone are 
sized using the proprietary MegaB ACE system, or an equivalent. Sizing is performed with 
multiple clones/lane. This multiplexing allows sizing to be performed in a cost and time efficient 
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manner. The multiplexing is performed with a liquid handling robot (e.g., Matrix PlateMate). 
After running the multiplexed fragments on MegaB ACE, and correlating the size of the fragment 
with the E. coli clone containing the insert, the fragments are analyzed to determine suitability 
for sequencing. 

Example 8. Comparison of clone complexity with and without use of a sizing step 

The effect of using a clone sizing step on the complexity, i.e., the representation of rarely 
transcripts, of the resulting clones, is shown in FIGS. 4A and 4B. In FIG. 4A, no sizing step was 
used, while clone sizing was used in the identification of the clones shown in FIG. 4B. Shown in 
the figures is a comparison of the frequencies (expressed in percentage) of clones derived from 
transcripts present at varying levels. The outer numbers represent the prevalence of a particular 
clone sequenced, and the inner numbers represents the percentages of the total number of clones 
sequenced that fall into this abundance class. As illustrated in FIG. 4A, the sequencing results 
that were obtained without the use of the sizing filter demonstrated that only a small percentage 
of the total number of fragments that were sequenced were included low copy number fragments 
(i.e., singletons, duplicates, and triplicates). Specifically, singletons were found to comprise only 
2% of the total number of fragments sequenced, while fragments that were present at greater than 
51 copies comprised 38% of the total fragments sequenced. In contrast, as illustrated in FIG. 4B, 
the sequencing results that were obtained with the use of the sizing filter were enriched for clones 
from low abundance transcripts (i.e., singletons, duplicates, and triplicates). These clones 
constituted approximately 33% of the total fragments sequenced. In contrast, without the use of 
this sizing filter, these fragments were found to only comprised a total of 8% of the sequencing 
results. 

Equivalents 

Although particular embodiments have been disclosed herein in detail, this has been done 
by way of example for purposes of illustration only, and is not intended to be limiting with 
respect to the scope of the appended claims that follow. In particular, it is contemplated by the 
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inventor that various substitutions, alterations, and modifications may be made to the invention 
without departing from the spirit and scope of the invention as defined by the claims. For 
example, the selection of the specific tissue(s) or cell line(s) that is to be utilized in the practice 
of the present invention is believed to be a matter of routine for a person of ordinary skill in the 
art with knowledge of the embodiments described herein. 
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WHAT IS CLAIMED IS: 

1 . A method of screening a population of nucleic acids for a novel sequence, the method 
comprising: 

providing a population of nucleic acid sequences; 

partitioning said population into one or more subpopulations of nucleic acids; 
identifying a first nucleic acid sequence in the subpopulation of nucleic acid sequences; 

and 

comparing the first nucleic acid sequence to a reference nucleic acid sequence or 
sequences, wherein the absence of the first nucleic acid sequence in the reference nucleic acid or 
nucleic acid sequences indicates the first nucleic acid is a novel nucleic acid sequence. 

2. The method of claim 1, wherein said DNA population is a cDNA population derived 
from a population of RNA molecules. 

3. The method of claim 2, further comprising partitioning the RNA molecules. 

4. The method of claim 2, wherein said cDNA population is derived from the 5' ends of the 
RNA molecules. 

5. The method of claim 2, wherein said cDNA population is derived from the interior 
regions of the RNA molecules. 

6. The method of claim 2, wherein said cDNA population is derived from the 3' ends of the 
DNA molecules. 

7. The method of claim 2, wherein said partitioning step comprises hybridization of a probe 
nucleic acid sequence to the population of nucleic acids. 
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8. The method of claim 2, wherein said partitioning step comprises digesting the cDNA 
molecules with one or more restriction enzymes. 

9. The method of claim 8, further comprising ligating adapter oligonucleotides to the 
termini of the digested cDNA molecules. 

10. The method of claim 9, further comprising amplifying the ligation products. 

1 1 . The method of claim 8, further comprising separating the amplified products. 

12. The method of claim 11, wherein said separating is by gel electrophoresis. 

1 3 . The method of claim 1 1 , wherein the first nucleic acid sequence is identified by 
comparing the size of one or more digestion products produced by a member of the 
subpopulation of nucleic acids to the sizes of fragments generated by the same restriction enzyme 
or enzymes in said reference nucleic acid or nucleic acids. 

14. The method of claim 11, further comprising 
recovering one or more size-separated digestion products; 
reamplifying the recovered products; and 

separating the reamplified products. 

15. The method of claim 14, wherein said separating is by gel electrophoresis. 

16. The method of claim 15, wherein the first nucleic acid sequence is identified by 
comparing the size of one or more digestion products produced by a member of the 
subpopulation of nucleic acids to the sizes of fragments generated by the same restriction enzyme 
or enzymes in said reference nucleic acid or nucleic acids. 
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1 7. The method of claim 9, further comprising: 

inserting the ligated adapter oligonucleotide into a cloning vector to form a vector-insert; 
transforming the vector-insert into a suitable host; 

culturing transformed host under conditions allowing for replication of the vector-insert; 
recovering the vector-insert from said host; and 

digesting the vector-insert with one or more restriction enzymes, thereby releasing said 
insert; and 

comparing the size of the insert to sizes of fragments generated by the same restriction 
enzyme or enzymes in said reference nucleic acid or nucleic acids. 

1 8. The method of claim 1 , wherein comparing is by determining at least a portion of the 
nucleotide sequence of the first nucleic acid sequence and comparing the nucleotide sequence to 
the nucleotide sequence of one or more reference nucleic acids. 

1 9. The method of claim 1 , wherein comparing is by hybridizing the first nucleic acid 
sequence to one or more of the reference nucleic acid sequences. 

20. A method for equalizing the representation of nucleic acids in a population of nucleic 
acids, the method comprising: 

providing a population of nucleic acid sequences, wherein said population comprises a 
first nucleic acid and a second nucleic acid having a nucleic acid sequence distinct from the first 
nucleic acid, and wherein said first nucleic acid is present at a higher level in said population 
than said second population; 

partitioning said population into one or more subpopulations of nucleic acids; and 

comparing the levels of said first nucleic acid sequence to the levels of said second 
nucleic acid sequence in the subpopulation of nucleic acid sequences, wherein a lower level of 
the first nucleic acid sequence relative to the second nucleic acid sequence indicates the 
representation of said first and second nucleic acid sequences are normalized. 
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21 . A method for producing a population of nucleic acid molecules enriched for 5' regions of 
mRNA molecules, the method comprising: 

providing a population of RNA molecules, said population including RNA molecules 
having a 5' terminal Gppp cap structure and a 5' terminal phosphate group; 

contacting said population of RNA molecules with a phosphatase under conditions that 
result in removal of the 5' terminal phosphate group while leaving the 5' terminal Gppp cap 
structure intact; 

inactivating said phosphatase; 

contacting the population of RNA molecules with a pyrophosphatase under conditions 
that result in the removal of the 5' terminal Gppp and the formation of a 5' phosphate group; 

annealing an oligonucleotide in the presence of an RNA ligase to form a hybrid molecule; and 

forming a cDNA from said oligonucleotide. 

22. A method of identifying an RNA sequence in a sample comprising a plurality of RNA 
sequences, the method comprising: 

synthesizing cDNA copies of a plurality of RNA species to form a cDNA sample; 

determining the size of one or more of said cDNA molecules in said cDNA sample; 

comparing the size of said sample with the size of a reference nucleic acid; and 

thereby identifying the cDNA sequence. 

23. The method of claim 22, wherein said cDNA molecules are digested with one or more 
restriction enzymes prior to the determining step. 

24. The method of claim 23, further comprising ligating adapter oligonucleotides to the 
termini of the digested cDNA molecules prior to the determining step. 
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25. The method of claim 22, wherein said identifying step comprises comparing the size of 
one or more digestion products produced by one or more said cDNA molecules to a reference 
nucleic acid or nucleic acids. 

26. A method of identifying an RNA sequence in a population of RNA sequences, the 
method comprising: 

(a) removing 5' terminal pppG from RNAs in said population to form a population of 
RNAs having terminal 5' phosphate groups; 

(b) ligating a linker oligonucleotide to the terminal 5' phosphate groups of RNA 
molecules in said population of RNAs; 

(c) synthesizing complementary cDNA molecules from said population of RNA 
molecules to form a cDNA sample; 

(d) digesting said complementary cDNA molecules with at least one restriction enzyme; 

(e) ligating an adapter molecule to the digested cDNA molecules; 

(f) amplifying the molecules produced in step (e); 

(g) identifying the amplified molecules of step (f); and 

(h) comparing the amplified molecules to one or more reference nucleic acids. 
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METHOD OF IDENTIFYING NUCLEIC ACIDS 



Disclosed are methods for identifying nucleic acids in a sample of nucleic acids in which 
nucleic acids are initially present in unequal amounts. The methods include partitioning the 
starting population of nucleic acids to form one or more subpopulations, and then identifying 
nucleic acids that are present in different amounts in the partitioned nucleic acid sample as 
compared to the starting population. 
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